Phase diagram of diluted magnetic semiconductor quantum wells 
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The phase diagram of diluted magnetic semiconductor quantum weUs is investigated. The inter- 
action between the carriers in the hole gas can lead to first order ferromagnetic transitions, which 
remain abrupt in applied fields. These transitions can be induced by magnetic fields or, in double- 
layer systems by electric fields. We make a number of precise experimental predictions for observing 
these first order phase transitions. 



Semiconductor and ferromagnetic materials have com- 
plementary properties in information processing and 
storage technologies. Recent advances in the MBE 
growing techniques have made possible the fabrication of 
Mn-based diluted magnetic semiconductors (DMS) with 
a rather high ferromagnetic-paramagnetic critical tem- 
perature Tc. 1^ In semiconductors it is possible to mod- 
ulate spatially the density of carriers by changing the 
doping profiles but in DMS it is also suitable to vary the 
magnetic order of the carriers by changing the magnetic 
ion densities. The combination of these two possibilities 
opens a rich field of applications for these materials. 

The high critical temperature DMS's have a high con- 
centration, c, of Min?^ ions, randomly located. The itin- 
erant carriers in the Gai-xMnxAs systems are holes and 
their density, c*, is much smaller than the magnetic ion 
density. In doped semiconductor the spin S=5/2 M-n?^ 
ions feel a long range ferromagnetic interaction created 
by the coupling mediated by the itinerant spin polarized 
carriers. |^-^ This interaction has some resemblance with 
the one observed in the pyrochlores where a sim- 

ilar type of coupling between a narrow electronic band 
and magnetic ions is assumed to exist. In the latter case, 
the phase diagram is significantly different from that of 
a conventional itinerant ferromagnet, showing first order 
transitions and phase separation Q and/or the forma- 
tion of localized textures near the transition temperature 
These features lead to interesting transport prop- 
erties, like colossal magnetoresistance. We expect these 
effects to be amplified in low density two-dimensional 
(2D) doped semiconductors, where the carrier-carrier in- 
teraction can play a significant role, and also favors fer- 
romagnetism pl| . 

In this work, we investigate the nature of the phase di- 
agram of diluted magnetic semiconductor quantum wells, 
with emphasis on the existence of abrupt transitions for 
experimentally accessible parameters. We analyze quan- 
tum wells made of GaMnAs growth in the < 001 > di- 
rection and with thickness w. The confinement of the 
carriers in the z-direction can be obtained by sandwich- 
ing the GaMnAs system between non-magnetic GaAlAs 
semiconductors. 

The system is described by the following Hamiltonian, 



H = Hm+ Hhoies + J ^ S/ • s, (5 (r, - R/) 



(1) 



We now analyze the three terms contributing to Ti. 

i) Ti.M is the Hamiltonian of the localized spins interacting 
with an external magnetic field B. Direct interactions 
between the magnetic moments of the Mn ions are much 
smaller than the interaction with the carrier spins and 
therefore, unlike in the case of the pyrochlores, we neglect 
this coupling. 

ii) Ti.hoies is the part of the Hamiltonian which describe 
the itinerant holes. It is the sum of the kinetic energy 
of the holes and the hole-hole interaction energy. We 
treat the kinetic energy of the carriers in the framework 
of the envelope function approximation. The confine- 
ment of the carriers in the quantum well produces the 
quantization of their motion and the appearance of sub- 
bands which, for sufficiently narrow quantum wells, have 
a clear light-hole or heavy-hole character. We assume 
that the hole density is low enough so that only one of 
the subbands, heavy-like, is occupied. With this, the in 
plane motion of the holes can be approximated by a sin- 
gle parabolic band of effective mass m* . We describe the 
interaction between the electrical carriers with the local 
spin density approximation (LSDA) fl^ ]. Except at very 
low densities, the LSDA gives an accurate description of 
the electron gas confined in quantum wells |p^ , p^ |. Since 
the hole g-factor is much smaller than the Mn (^-factor 
we neglect in Tihoies the coupling between the hole spins 
and the applied magnetic field. 

iii) The last term is the antiferromagnetic exchange inter- 
action between the spin of the Mn^+ ions located at R/ 
and the spins , s^, of the itinerant carriers. The interac- 
tion between ions mediated by the conduction holes is of 
long range. Thus we will assume that thermal distribu- 
tion of the orientation of the Mn spins is that induced 
by an effective field, due to the holes, which should be 
calculated selfconsistently. 

In order to obtain the phase diagram of the DMS sys- 
tem we have to minimize the free energy T. The criti- 
cal temperatures for the ferromagnetic to paramagnetic 
transition in DMS is typically smaller than lOOif , and 
for these temperatures we can consider that the electron 
gas is degenerate. Hence, the only temperature depen- 
dence in J- is due to thermal fiuctuations of the Mn spins. 
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Treating the holes in the mean field approximation the 
free energy per unit area takes the form: 

T = T,ons + ^\nla{l + e) + EUn2D,0 ■ (2) 

Here ^ is the carrier spin polarization, n2D~c* w is the 
two dimensional density of holes in the system, E^c is the 
exchange correlation energy for the holes and J-ions is the 
contribution of the magnetic ions to the free energy: 
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the effective magnetic field than the M n spins feel. In ob- 
taining Eq.|^] we have assumed that the hole wave func- 
tion in the z-direction has the form w~^/'^ . 

The phase diagram with parameters J — 0.15eV nm'^, 
ion concentration c = Inm^'^ and width of the well, d = 
lOnm, is shown in fig.{|^}. We include a single hole band 
of effective mass m\\ = O.llme, and a dielectric constant 
eo = 12.2 [|. 



novelty of our calculation is the identification of a first 
order transition to a fully polarized state at intermedi- 
ate densities. This transition takes place at higher tem- 
peratures than that at which the magnetic susceptibility 
of the system diverges. The existence of a first order 
transition between the paramagnetic and the ferromag- 
netic phases implies the existence of a discontinuity in the 
chemical potential at Tc- This discontinuity indicates the 
existence of a region in the phase diagram where phase 
separation can take place. At higher densities, when a 
continuous transition between ^ = and £, ^ occurs, 
we also find a first order phase transition between the 
partially polarized system, ^ < 1, and the fully polarized 
phase, C=l- 

In order to understand better the phase diagram it is il- 
lustrative to derive the existence of these first order tran- 
sitions analytically. For doing that we treat the hole-hole 
interaction in the Hartree-Fock (HF) approximation. We 
obtain that the HF results are in rather good agreement 
with the obtained using the LSDA. In the HF approach 
the hole energy per unit area can be written as, [[L8[ 
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FIG. 1. Phase diagram of a magnetic semiconductor quan- 
tum well, using the parameters described in the text. Full 
line: first order transition. Broken line: second order phase 
transition. The holes in that phase to the left of the full line 
are completely polarized. The inset shows the polarization of 
the holes for kp ~ 0.2nm~^ (broken line) and kp = 0.39nm~^. 

The calculations have been done by minimizing the 
value of the free energy, Eq. |^ and using the expression 
given by Vosko et al. for the exchange correlation po- 
tential. In the high density region the line of continu- 
ous transitions agrees with that obtained from the di- 
vergences of the magnetic susceptibility H. The main 



where the last term is the exchange energy of the holes. 
The free energy of the ions is given by Eq.|^. Near a 
paramagnetic-ferromagnetic transition, the hole spin po- 
larization is small and we can expand the total free energy 
in terms of powers of ^: 
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with 



We expect a paramagnetic-ferromagnetic transition when 
the quadratic term in ^ is zero. This implies a critical 
temperature 



Tc = 
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At high densities, this approximation gives Tc ~ 8.7K, 
in good agreement with the LSDA calculation, and with 
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the RKKY solution Note that the only dependence 
of Tc on n2d is through the hole-hole interaction. This is 
due to the fact that the two-dimensional density of states 
is energy independent. In three dimensional systems the 
kinetic energy scales as n^/'^ and Tc is proportional to 
v}/^ in agreement with the RKKY theory. The order of 
the transition can be inferred from the sign of the quartic 
term in Eq. If the quartic term is positive the tran- 
sition is of second order, while a negative quartic term 
implies the existence of a first order phase transition. Us- 
ing the previous expression for Tq, a first order transition 
takes place if: 
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which can be solved to give n^^J'^^ « 2.1 x lO^^cm^^. 
This estimate is also in good agreement with the LSDA 
results shown in fig.{|l|}. 

Some comments about the LSDA exchange correlation 
potential are in order: 

i) In the HF treatment, the existence of a negative quartic 
term in the expansion of the exchange energy in powers 
of ^ is the source for the appearance of first order phase 
transitions. In the LSDA, the intermediate spin polar- 
ization correlation energy is obtained [ p^ , p^p^Jl5[ | by as- 
suming that it has the same polarization dependence as 
the exchange energy, so that the LSDA expression for 
Exc also leads to a negative quartic term when expanded 
in powers of ^, and a first order transition is expected. 
On the other hand, correlation effects ||l^ weaken the 
spin dependence of the interaction energy as compared 
with the results in HF, and we find that the LSDA gives 
a value for n^^J''^ slightly smaller than the observed in 
the HF approximation. Numerical evaluation [Q of the 
partially spin polarized Exc also shows a negative quartic 
term in the expansion of Exc versus ^. Hence, we believe 
that the existence of a first order phase transition is a ro- 
bust result, independent of the model used for describing 
the interaction between carriers. Note, finally, that the 
existence of a first order transition implies the absence of 
long range critical fluctuations, lending further support 
to the adequacy of the methods used in this work. 

ii) The interpolation formula used for describing the ^ de- 
pendence of Exc is not analytic at ^^=1, Eq.Q. This 
implies that there is, for each pair of values n2d and T, a 
^ range near ^=1 which cannot be reached by minimiz- 
ing the total free energy. When dT"^ jd^ is smaller than 
zero, the system prefers to be completely spin polarized. 
More accurate descriptions of the LSDA Exc [0 results 
in an analytic behavior at ^=1. Therefore we believe 
that, in the results shown in fig.{|l|}, the discontinuous 
transition which occurs at lower temperatures than the 



second order transition T^, is, probably, a spurious result 
due to the use of a HF-like interpolation for Exc- 
We remark again here that the appearance of a first order 
transition when decreasing n2d is a real result, not related 
with the non analyticity of the LSDA expression for Exc- 

Now we analyze the effect in the phase diagram of an 
external magnetic field. For B ^^0, the discontinuous 
transitions shown in fig.{|l|} are changed into transitions 
between phases with two different polarizations. General 
thermodynamic arguments show that this line of first or- 
der transitions should end in a critical point, in an anal- 
ogous way to the liquid-vapor phase diagram. This crit- 
ical point belongs to the Ising universality class. In the 
present calculations, the first order transition persists at 
all fields. This is due to the non analytic of the exchange 
energy at ^=1. We expect that a correct description of 
Exc would induce the termination of the line of first order 
transitions at a critical value of the field. Finite temper- 
ature effects in the hole gas, not taken into account, also 
will weaken the ^ dependence of the free energy, leading 
to the existence of a critical point. 

The line of discontinuous transitions in an applied field 
is shown in fig.{0}. 
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FIG. 2. Line of first order phase transitions as function of 
applied field and temperature, for fixed density. 

The existence of discontinuous transitions leads to the 
possibility of phase separation. Using Maxwell construc- 
tion, we find a region of phase separation near the line of 
first order transitions. In order to analyze its occurrence, 
we need to include electrostatic effects. The simplest sit- 
uation where phase separation can be observed is in a 
bilayer. Let us imagine two wells with nominal hole den- 
sity equal to riQ. At the temperatures and fields where the 
first order transition occurs, there are two phases with the 
same free energies, Txin^) = Ti(n^- The chemical po- 
tential of these two phases differ at this point, and we de- 
fine A/i = /ii — /i2- This difference in chemical potentials 
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will induce a transfer of charge between the two layers, 
5n, until A/i is compensated by the induced electrostatic 
potential, V « {e'^6nd)/eo, where d is the distance be- 
tween the layers. Thus, we obtain Sn « (A/ieo)/(e^rf). 
For reasonable values of c? 10 — 50nm, we find that the 
charge transfer is small, Sn/no ^ 10^^. The change in 
the electrostatic barrier, V « A/j, is, however, compara- 
ble to the width of the band, and can change significantly 
the transport properties. 

In a bilayer system, the first order transitions analyzed 
here can be induced by an applied electric field. The field 
induces a difference in the chemical potential of the two 
layers, which leads to a charge transfer. By suitably tun- 
ing the parameters, the density in one of the layers will 
reach the value at which the first order transition dis- 
cussed above takes place. At this point, there will be an 
abrupt change in the charge transfer, which can be mea- 
sured with standard capacitive techniques ||l6|,|r^. The 
charge transfer as function of electric field, for reasonable 
parameters, is shown in fig.{^. 



LB thanks A.H.MacDonald, C.Tejedor and A.Rubio 
for helpful discussions. We acknowledge financial support 
from grants PB96-0875 and PB96-0085 (MEC, Spain) 
and (07N/0045/98 and 07N/0027/99) (C. Madrid). 




2 3 

V+ (meV) 

FIG. 3. Charge transfer in a bilayer system as function of 
the applied electric field. 

In conclusion, we have analyzed the possible dis- 
continuous transitions in two-dimensional diluted mag- 
netic semiconductors. We find that the interaction be- 
tween carriers can lead to first order transitions in two- 
dimensional quantum wells. At the transition, the holes 
become fully polarized. This transition can be induced 
by a change in the density of the hole gas, the temper- 
ature, a magnetic field, and, in multilayer systems, an 
applied electric field. At these transitions, the minority 
spin band becomes empty, and the density of states at the 
Fermi level is reduced by one half. In double layer sys- 
tems, significant electrostatic barriers can appear near 
the transition. This change can alter significantly the 
transport properties pO{ . Thus, it can be important in 
the operation of devices made with these materials. 
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